— - -fSTECTS c:* li-AVV IONS on ELECTRON ^*«WKF-\TU5L5 


IN THE SO’ 

C* ^ ~ 


Curonj. and solar wind 


NAS.V-Gcddard Space Flight Center 
Greenbelt- Maryland 20771 


ABSTRACT 
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The effects of the reduction in the thermal 
conductivity due to heavy ions on electron temperatures 
in the solar corona and solar wind are examined. Large 
enhancements of heavy ions in the corona appear to be 
necessary to give appreciable changes in the thermal 
gradient of the electrons. These enhancements, if they 
should occur, may contribute to the understanding of 
some low values of solar wind temperature measurements 
at 1 AU. 
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I. INTRODUCTION 


At 1 AU, observations indicate that solar wind ions 
have nearly the same mean flow velocities and that their 
elemental abundances are very approximately the same as that 
of the photosphere (Formisano and Moreno , 1971) . The 
abundances are such that the consideration of heavy ions 
may be expected to give thermal conductivities that are 
almost the same as those obtained for only protons and 
electrons. In their review, Formisano and Moreno find 
that many diffusive models give mean flow velocities for 
the heavy ions that are less than mean flow velocities of 
H + near the sun. If these differences in flow velocities 
and thus the enhancement of heavy ions in the corona should 
occur, thermal conductivities of the electrons may be modified 
enough to give appreciable changes in the temperature 
gradient near the sun. The purpose of this study is to 
examine the effects of such possible enhancements on 
electron temperatures in the corona and interplanetary 
space . 

Many simplifying assumptions are made. Spherical 
symmetry and time independence are assumed. Heat sources 
are assumed cnly at the base of the corona. Collisional 
heat exchange between ions and electrons is usually neglected. 
Magnetic field and viscosity effects are not included. 

In the Tollowing section, the electron equations 
and parameters are discussed. Results are given and 
discussed in Section III. 
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II. EQUATIONS 

The electron energy equation that is appropriate for 
the assumed conditions is: 

2k — -— — — ( r 2 KT 5 ^ 2 — ^ = CL' (1) 

2 dr n dr nvr 2 dr \ dr/ 

where T, n, and v are the electron temperature, density, 
and mean flow velocity. CL 1 is a term that gives collisional 
heat transfer, with ions. The same equation is obtained 
for distribution functions expressed in terms of the mean 
mass flow velocity or mean flow velocity of the electrons 
( Burgers , 1969) . The procedures for obtaining K are, 
however, quite different. 

K values used here are obtained from an adaptation 
of the derivations by Burgers and are compared with values 
obtained by Spitzer and Harm (1953) . Consider first the 
derivation where equations are given in terms of the mean 
mass flow velocities. These equations appear to be applicable, 
although differences in flow velocities of various ions 
relative to the mean mass flow velocity may be quite large, 
since these differences appear to be small compared with mean 
thermal velocities of electrons except for T«10 I *K. 

In his treatment of heat flow. Burgers uses the momentum, 
heat flow, and electric current equations for a plasma with a 
single ion species with charge, Z, but omits the pressure qradient. 

For velocities relative to the mean mass flow velocity, heat 
transfer by ions is considered small compared to that by 
electrons. With no net current flow, the following K may 
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be obtained from his results; 



( 2 ) 


The inclusion of the pressure gradient appears to give the 
same result. For the same assumptions, this result may 
be extended to mixtures of ions by replacing Z with 
^ n i Z i/^ n i Z i w ^ ere the summations are over the ions. 

For equations based on the mean flow velocity of the 
electrons, the derivation is somewhat simpler. The assumptions 
of little heat transport by ions for velocities relative 
to the mean mass flow velocity, and no net current flow 
may be used with Burgers* Equations (2.17a) and (32.13) to 
give the same K as above. 

Spitzer and Harm appear to have used methods equivalent 
to the ones based on the mean mass flow velocity. Their 
value of K (a similar K is used by Hartle and Sturrock , 1968) 
is about 2.4 times larger for Z=l. A part of the difference 
may be due to the smaller thermoelectric field that Spitzer 
and Harm obtained compared to that obtained by Chapman (1958) 
and Burgers. Decreases in K due to increases in Z are 
somewhat greater for Burgers' K than for Spitzer and Harm's. 

Equation (1) may be transformed by using: 


d&ogT m 
dZogr 


(3) 


d&ogn _ 
dJtogr 


e 


(4) 


to give: 
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aiffF - + < { -i e+CL > - tkHf 

CL gives heat transfer due to collision with ions and for 
T in 10 ®K, R in solar radii, and average conditions in 
the solar wind is approximately: 


CL = -E 


T - T 


0.9 


5 1 2 


(5) 


RT 3/ * (vxl0~ 7 ) 2 

CL is not included in the calculations below; some effects 
due to its inclusion are discussed. 

For the solar wind electron flux in Equation (la) , a 
median value of 2.6x10 8 per cm 2 sec at 1 ?U is used (Hundhausen 
et al . , 1970). Various values of Z[R) are used and are 
discussed below. e(R) that is compatible with electron 
density measurements is used. From the continuity equation 
e 5s: 


e = 2 + 


dtnv 

<PnR 


( 6 ) 


For 

d&nv ^ 7.0 

dJlnR R xn 

n values of 3.2x10* per cm 3 at R=1 , 6.5 per cm 3 at 1 AU and 
values at other R that are near measured values are possible. 
T results appear to be quite insensitive to the that is 

assumed provided that it gives n that are near coronal and 
interplanetary densities. 
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Values of the electric field are also approximately 
evaluated. Leaving out the small terms in the electron 
momentum equation, the electric field, E, is given by: 

dJln(nT) eE 
dr ” kT 


eE may be expressed in terms of the gravitational force: 


eE 


m*m 

P 


g 


where m* is in units of the proton mass. With T in 10 6 K 
and R in solar radii, m* is: 

★ TP 

"■ =20 (5+E) (9) 

Near the sun where flow velocities are small, hydrostatic 
equilibrium is approached; m* is then approximately the 
average particle mass including electrons, m* values near 
the sun should be 0.5 for a H + and e c r rona and should be 
>0.5 when other heavy ions are present. 
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III. RESULTS 

Once Z (R) is specified, simultaneous solutions to 
Equations (la) and (3) are most readily obtained by specifying 
T and 6 at some large R and integrating to R=l. For a 
given 6 at the outer boundary, T at the boundary is obtained 
from Equation (la) with the two derivative terms, which may 
be expected to be quite small, set equal to zero. 

A number of solutions, obtained through a computer 
for Z=1 (two fluid model), are shown in Figure 1. The 
various solutions are labeled with 6 at the outer boundary. 

The x symbols give the range of T that have been measured 
by Montgomery et al . (1968); the circles give the range 
of T that Serbu (1972) has measured. 

The curve with initial 6=0.30 is quite near results 
obtained by Hartle and Sturrock (1968) . This T result and 
their Tp result may be expected to give near minimum solar 
wind velocities. Solutions with higher T than this one 
may be expected to give solar wind velocities with T 
values near those of Hartle and Sturrock but tend to give 
rather high values of T near 1 AU. Solutions with lower 
T probably require Tp values that are higher than those 
of Hartle and Sturrock in order to give solar wind velocities. 

Solution above the 6=0.4 one are concave upward and 
may be expected to approach 6=2/7 at very large R. Solutions 
below the 6=0.4 solution are concave downward for R>10. If 
the collision term is not considered, 6 is expected to increase 
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with R for these solutions and T can become zero at some 
finite R. As T becomes small, however, the collision 
term (Equation (5)) can become appreciable; T solutions 
that approach at large R appear to be possible. At 
these small T and n values, other considerations such as 
viscosity and deviations from spherical symmetry may be 
the important ones in determining T. 

Solutions for other solar wind fluxes may be obtained 
from those of Figure 1 by moving solutions up or down. 

This is possible since the last term of Equation (la) and 
effects due to variations in JlnA are small. Solutions with 
increases (decreases) in fluxes by some factor are approxi- 

2/5 

mately given by shifting the curves up (down) by the (factor) ' . 

Solutions with a factor of 2 increase in flux are approximately 
given by moving all curves to higher T by 2 2 ^ 5 or 1.32. 

Similarly solutions for a constant composition with K 
decreased because of Z>1 may be obtained by a similar 
translation. For Z=1.21, which is near the value for the 
average solar wind at 1 AU the solutions are approximately 
obtained by shifting the solutions of Figure 1 upward by 
(2 JS/5 + 1.3x1.21) 2/ 7(2v^?5 + 1.3x1) 2/5 or 1.05. This is 
such a small amount that it appears to be quite certain 
that solutions for the average composition of the solar wind 
at 1 AU with no enhancement should be very near those for 
two fluid models. 
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The changes in the solution owing to enhancements of 
heavy ions near the sun are investigated by using the 
following : 


Z(R) = 1.21 + - 2 — - ~ 1,21 — ( 10 ) 

1+ (JtogR/i,og4) 3 

Solutions for 6=0.5 at R=10 4 and Z (1) values of 1.21 (lowest 
curve), 2, 4, and 8 (highest curve) are given in Figure 2. 

Other solutions with different 6 at the outer boundary 
give similar curves; increases in T at R=1 are by almost the 
same factors. 

It appears that reductions in the thermal conductivities 
due to the enhancement of heavy ions in the corona can 
contribute to the reduction of calculated electron temperatures 
at 1 AU while permitting high temperatures near the sun. 

The magnitude of these effects can be much less than or 
much greater than is shown in Figure 2 depending on the R 
dependence of Z. 

The amount of enhancements of heavy ions that are necessary to 

qive Z(l) values are illustrated in Table 1. The 0 in column 1 represents 
elements from C through Mg; S is for all elements heavier 
than Al. Column 2 gives assumed average charges. Column 3 
gives assumed relative densities at 1 AU. Columns 4-7 
give assumed enhancements at R=l. Below each of the last 
4 columns are Z(l) and hydrostatic m** values for each 
set of assumed enhancements. These assumed enhancements 
are large; however, enhancements derived from mean flow 
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velocities of heavy ions relative to hydrogen have been 
found to be quite Large for diffusive coronal models 
( Geiss et al, 1970, Nakada, 1970). The conditions assumed 
in their models, however, may be quite different from the 
real corona. In a different diffusive model, Alloucherie 
(1970) does not find appreciably different mean flow 
velocities near the sun. 

The m** in Table 1 is the approximate hydrostatic 
value and depends only on the composition at the base of 
the corona. It is approximately: 


In. Z .m. 

m** = — 1 x . x - (11) 

In. Z. 2 
i i 

where the summations are over the ions and electrons and m^ 


are in units of m^. This m** is compared below with the m* 
(Equation (9)) that is obtained from the electron distribution. 


In Figure 1, the solution with the 6=0.30 label and 
which is near results of Hartle and Sturrock has m*(l)=1.02. 


This could indicate th^ presence of a considerable amount of 
heavy ions for the electron distribution that was assumed, 
m* is, however, very sensitive to the e or the electron 
distribution that is assumed. By keeping the density at 1 AU 
the same and reducing the density at the base of the corona, 
it is possible to obtain m*(l)=0.5 and thus consistency with 
the two fluid models with only minor changes in the T(R) 
solution. 


m* values for the solutions of Figure 2 are shown in 
Figure 3. Only near the sun where ion velocities may be 
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expected to be small is there a meaningful comparison with 
the m**(l) of Table 1. Because of the sensitivity of the 
m* results to the electron distribution, the rather weak 
conclusion that can be drawn from this comparison is that 
it appears possible to obtain consistency ' >* cwet n** and m* 
with considerable enhancements of heavy ions and still have 
an electron distribution that is near measured values. 

Some trends in proton temperatures near the sun due 
to the presence of heavy ions may be obtained from the 
proton energy equation. Neglecting certain small terms 
Burgers' Equation (33.8) may be written: 



= m^g (**l+m*+W) 


(12) 


W is a term that gives heat transfer between various species 
and protons; for this discussion, the main contribution to 
W is considered to be due to differences in mean flow 
velocities of protons and heavy ions and not due to temperature 
differences between the ions. 

For the two fluid models, W is zero and the electric 
field is given by m^=o.5. The main reason for the drop 
in Tp near the sun is work done against the gravitational 
field minus the work done on the protons by E. If W were 
to remain zero and m* increase, a slower drop in T p compared 
to the two fluid models may be expected. Consider, for 
example, the situation with a near average solar wind 
composition and no enhancement in the corona (column 4 of 
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Table 1.). The flow velocities of all ions would be the 
same so W-0. But with m* =0.70, the net work done per unit 
distance wculd be reduced by 40% so that T may be expected 

ir 

to drop more slowly with radial distance than it does for two 
fluid models. If m*(l) were as large as 1.0 with small W, 

Tp may be expected to remain almost constant with distance 
near the sun. 

For enhancements that give m*(l) ,.ear 1.0, it is likely 
that W should be appreciable and negative. Two assumptions 
are made to illustrate the effects of appreciable W. The 
proton momentum equation is used to evaluate a term that is 
almost W with the assumption that np/n e is constant with r. 
This term is substituted into the proton energy equation with 
the assumption that the flow velocity of H is large compared 
with flow velocities of He ++ and other heavy ions. This gives 
the result that Tp can be constant with radial . .ance near 
the sun for m* between 1.1 and 1.4. This range is due to 
assumptions about the predominant heavy ion. Enhancements 
of heavy ions thus appear to be able to give Tp near the sun 
that are considerably highor than is given by the two fluid 
models. Tp greater than electron T for some distance near 
the sun appear to be possible. For sufficiently large m*, 

Tp may increase with distance near the sun. 

The results of this study are that enhancements of 
heavy ions in the corona tend to give lower electron 
temperatures and higher proton temperatures at 1 AU when 
compared with Hartle and Sturrock's two fluid models of 
the sclar wind. Considering the many simplifications and 
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that the momentum and energy equations for all of the ions 
are not yet solved and may even be insoluable, it can only 
be concluded that enhancements may contribute to the under- 
standing of solar wind temperature measurements but are n* ' 
necessarily required. 





TABLE 1 


Average Charges, Z(l), and Electric Fields Expressed in 
Terms of Near Average Particle Masses, m**(l), for Assumed 
Enhancements Near the Sun. 


Element Charge Relative Densities Assumed Enhancements 

at 1 AU 


H 

1 

1.0 

1 

1 

1 

1 

He 

2 

.05 

1 

10 

20 

40 

0 

8 

.002 

1 

10 

100 

400 

S 

11 

.0002 

1 

10 

100 

1000 



2(1) 

1.21 

2.07 

4.2 

6.2 



m**(l) 

0.70 

1.07 

1.66 

1.94 
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FIGURE CAPTIONS 


Figure 1. Solutions to the electron energy equation 
for two fluids for various values of 
6=-d£nT/di,nR at the outer boundary. 

Figure 2. Solutions to the electron energy equation with 

the thermal conductivity modified by enhancement 
of heavy ions near the sun. Assumed values of 
average charge, Z(l), at the base of corona are 
1.21 (lowest curve), 2, 4, and 8 (highest curve). 

Figure 3. Electric fields that are obtained from the electron 
density distribution and the solutions of Figure 2. 
m* gives the ratio of the electric field force 
to the magnitude of the gravitation force on 
protons . 







